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We study the thermodynamical properties of electrically charged black hole solutions of a nonlinear 
electrodynamics theory defined by a power p of the Maxwell invariant, which is coupled to Einstein 
gravity in four and higher spacetime dimensions. Depending on the range of the parameter p, 
these solutions present different asymptotic behaviors. We compute the Euclidean action with the 
appropriate boundary term in the grand canonical ensemble. The thermodynamical quantities are 
identified and in particular, the mass and the charge are shown to be finite for all classes of solutions. 
Interestingly, a generalized Smarr formula is derived and it is shown that this latter encodes perfectly 
the different asymptotic behaviors of the black hole solutions. The local stability is analyzed by 
computing the heat capacity and the electrical permittivity and we find that a set of small black 
holes are locally stable. In contrast to the standard Reissner-Nordstrom solution, there is a first- 
order phase transition between a class of these non-linear charged black holes and the Minkowski 
spacetime. 



I. INTRODUCTION 

The nonlinear electrodynamics models have been proved to be excellent laboratories in order to circumvent some 
problems that occur in the standard Maxwell theory. Indeed, the interest for nonlinear electrodynamics has started 
with the precursor work of Born and Infeld, whose main motivation was to modify the standard Maxwell theory in 
order to eliminate the problem of infinite energy of the electron [IJ. However, due to the fact that the Born-Infeld 
model did not fulfill all the hope, nonlinear electrodynamics theories become less popular. The recent renewal interest 
on the nonlinear electrodynamics theories is essentially due to their emergence in the context of low-energy limit of 
heterotic string theory, where a Gauss-Bonnet term coupled to quartic contractions of Maxwell field strength appear, 
and black hole solutions can be obtained It is also important to mention that nonlinear electrodynamics theories 
are a powerful tool for the construction of regular black hole solutions [3] . 

The thermodynamics properties of nonlinear electrodynamics theories is also an active research area in the current 
literature, see e.g. [4] in the case of the Born-Infcld model. A very appealing property which is common to all 
the nonlinear electrodynamic models lie in the fact that these theories satisfy the zeroth and first law of black hole 
mechanics [S]. This property renders more attractive the studies of nonlinear electrodynamics models. However, in 
contrast with the standard Maxwell theory, the formula of the total mass, the so-called Smarr formula P , does not 
hold a priori for nonlinear electrodynamic theories [5|. Indeed, in the case of the Einstein-Maxwell system, the Smarr 
formula and the first law of black holes mechanics are closely related and each of them can be derived from the other. 
This correspondence between both formulas is due to the homogeneity of the mass in terms of the area of the horizon 
and the electric charge |7l[8]. However, in the case of nonlinear electrodynamic theories, this homogeneity property is 
not longer present and there is a priori no good reason to obtain a Smarr formula. Some attempts to generalize the 
Smarr formula in order to be in accordance with the variations expressed by the first law of black holes mechanics have 
shown to be possible for some particular magnetic solutions in four dimensions [S]. For completeness, we stress that 
the derivation of the Smarr formula as well as of the first law of thermodynamics are usually performed successfully for 
asymptotically flat spacetime through the Komar integrals [lOj . However, their extension to rotating asymptotically 
anti-de Sitter black hole solutions is not a simple task as it can be seen by reading the current literature on the topic 
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In the present paper, we consider higher-dimensional gravity coupled to a nonlinear electrodynamic source given 
by a rational power p of the Maxwell invariant. This theory exhibits electrically charged black hole solutions as 
long as the exponent p belongs to the set of rational number with odd denominator |13l 114] . Similar considerations 
and motivations were previously analyzed in three dimensions in |15j . Recently, this nonlinear model has attracted 
attention, and it has been used for obtaining Lovelock black holes [16], Ricci flat rotating black branes :17], magnetic 
strings [TH^ and topological black holes in Gauss-Bonnet gravity [T3] . A generalization of this model for the non abelian 
case was considered in [20J. The black hole solutions derived in |14j are divided in four ranges for the exponent p, and 
each of them correspond to a different spacelike asymptotic behavior of the metric. Among these different classes, 
there exist black hole solutions with unusual asymptotic properties as those that go asymptotically to the Minkowski 
spacetime slower than the Schwarzschild spacetime. There also exists a range for p for which the solutions are not 
asymptotically flat and, their asymptotic behavior is shown to grow slower that the Schwarzschild de Sitter spacetime. 
Extremal black hole solutions generalizing the extremal Reissner-Nordstrom solution arc also known for this model 



The main objective here is to provide a detailed analysis of the thermodynamical properties of these black hole 
solutions. In particular, we would like to explore the relation between the thermodynamical properties and the 
different asymptotic behaviors of the solutions. The thermodynamics of the electrically charged black hole solutions 
will be performed through the Regge-Teitelboim approach j21l . in which an appropriate boundary term is added to 
the Euclidean action such that the total action presents an extremum. On the other hand, owing to the fact that 
the Euclidean action is related to the Gibbs free energy, the identification of the mass and the charge, as well as the 
study of the global stability of the black holes will be considerably facilitated. We will derive a generalized Smarr 
formula by using the explicit expressions of the temperature and the entropy. Interesting enough, we will show that 
the different asymptotic behaviors of the black hole solutions are reflected through this Smarr formula. In addition, we 
will present two other different ways of obtaining of the Smarr formula. The first derivation is operated through the 
Komar integrals for all the class of solutions even those that are not asymptotically fiat. The same formula will also be 
obtained with the use of a Noether conserved current which is associated to a scale symmetry of the reduced action. 
The local stability of the black holes are analyzed by computing the heat capacity and the electrical permittivity. For 
a certain range of the exponent p, small black hole solutions with positive and negative mass will be shown to be 
locally stable. Finally, the global stability is studied through the Gibbs free energy in order to determine whether the 
electrically charged black hole solutions are most likely than the Minkowski background. As a result, we will establish 
that there exists a phase transition for a certain range of the exponent p. 

The plan of the paper is organized as follows. In the next section, we review the nonlinear electrodynamics model 
and its general black hole solutions [14]. In Sec. Ill, the thermodynamics properties of the system are studied through 
the Euclidean Hamiltonian formalism. The mass M and the charge Q are explicitly identified and a generalized Smarr 
formula is derived. In Sec. IV the local and global stability of these black holes is analyzed, while the last section is 
devoted to the conclusions and further prospects. Finally, two appendices are devoted to the different derivations of 
the Smarr formula. 



In |13l I14j . a nonlinear electrodynamics coupled to gravity in d spacetime dimensions^ was considered. This model 
is described by the action 



where k denotes the gravitational constant, a the coupling constant for the electrodynamical action, F^i, — dfj^A^, — 
d^A^ is the field strength and p is a rational number whose range will be fixed later. The field equations obtained by 
varying this action are given by 
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ll. NONLINEAR ELECTRODYNAMICS AND BLACK HOLE SOLUTIONS 




(1) 



G^, = 4/ta pF^pF/ F'P-^ - -g^, F^ 



(2a) 




(2b) 



Along the present article we consider d > 3. 
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Type 


b 


P 


Remarks 


I 


b>d-3 


1/2 <p< {d-l)/2 


Standard asymptotically flat case 


II 


< 6 < d- 3 


p> (d-l)/2 orp<-l/(d-4) 


Electric term with relaxed fall-off 


III 


6 < 


-l/(d-4) <p< 


Asymptotically non-flat case 


Log 


b = d-3 


p = (d- l)/2 with d odd 


Logarithmic case 



TABLE I: Classification of the solutions depending on its asymptotic behavior. The cases are labeled in the first column. In 
the second column the different cases are defined by the parameter b and, equivalently by p, in the third column. In the last 
column we describe the main spacelike asymptotic feature for each case. 



where F — Fa/sF"^ is the Maxwell invariant. The most general spherically symmetric solution with a radial electric 
field was found in [13] for the conformal case p = d/A, and in jl4] for p belonging to the set of rational numbers with 
odd denominators^. The general solution is described by the line element 

ds^ = -N\r)fir)dt' +4^+ r'dnl_,, (3) 



Pir) 



with 





= 1, 


fir) 


= 1 

r 




C 


Ftr — 


d-2 ' 
'p2p — l 



(4a) 

d^ + -b' (4b) 



(4c) 



where for convenience we have defined 



2Ka{~l)PC^P2P{2p-l)^ 2{pd-4p+l) 

{d-2){d-2p-l) ' 2p-l ' ^' 

These solutions contain black hole configurations depending on the values of d and p, and the integration constants 
A and C. 

Apart from the restriction of being a rational numbers with odd denominators, the exponent p can not belong to the 
set (0, 1/2] since in this case the scalar curvature of the solutions diverges at the infinity. The form of the lapse function 



suggests a natural partition depending on the exponent h that appears in ( 4b ) . These different ranges are given by 
6>d— 3, 0<6<(i — 3,6<0, and h — d—'i, which correspond to the different spacelike asymptotic behaviors of the 
metric, as it is shown in Table I. Solutions of type I have a behavior similar to the standard Reissner-Nordstrom one. 
By similar, we mean that the charge term in the metric decays faster than the mass term in the asymptotic region. 
At the opposite, solutions of type II correspond to black hole solutions which go asymptotically to the Minkowski 
spacetime slower than the Schwarzschild spacetime. In the case III, the black hole solutions are not asymptotically 
flat, and their asymptotic behavior is shown to grow slower that the Schwarzschild de Sitter spacetime. Finally, the 
critical value h — d — 3, which can only occur for odd dimensions d — 2p -\- 1, corresponds to the limit between the 
solutions which resemble to the standard Reissner-Nordstrom solution (type I) and those for which the asymptotic 
decaying to Minkowski spacetime are slower than the Schwarzschild spacetime (type II). The spherically symmetric 
solution in this special case involves a logarithmic dependance on the radial coordinate and is given by 

Ftr = (6a) 
r 

nr) = i--^,+^a{~m'^'c-'^,- (6b) 



^ This restriction on p arises from considering spherically symmetric real solutions with a purely electric radial field. 



4 



III. THERMODYNAMICS 



It is well known that the partition function for a thermodynamical ensemble can be identified with the Euclidean 
path integral in the saddle point approximation around the Euclidean continuation of the classical solution |22j . In this 
case, the Euclidean action Ie evaluated on the classical solution is related to the free energy G of a thermodynamical 
ensemble by Ie = PG, where f3 is the inverse of the temperature which corresponds to the period of the Euclidean 
time T. 

As a first step we write the action in the Hamiltonian form, and since we are concerned only with the static, 
spherically symmetric case without magnetic field, it is enough to consider a reduced action principle. The class of 
the Euclidean metric and electric potential to be considered are given by 



dr- 

JW 



r^dn 



2 

d-2) 



A = A{r)dT, 



where the radial coordinate r belongs to [r^,oo). In this case, the Euclidean reduced action obtained from ([ij reads 



Ie = -P nd^2 



2p 
2p-l 



d-2 



{2p- l)aiV(-2)55fer / V \ 



-Nr' 



d-2 



{PY 



3 



(1-/^) 



+ ^V'\+ K, (7) 



where V = AapNoo r'^^'^FP~^F^^ is the rescaled canonical radial momentum, (j) = A{r) is the electrostatic potential, 
is the area of the d— 2-dimensional unit sphere. The boundary term K appearing in ^ will be fixed by requiring 
that the action has an extremum on-shell [21 . Moreover, owing to the fact that the Hamiltonian action is a linear 
combination of the constraints, the value on the action on-shell is given by the boundary term K. 

The equations of motion obtained by varying the reduced action with respect to N , f'^,V and (j) are given by 



d-3 



r r 
N' = 0, 



, (i-,n = 



2Ka{2p- 1) (-2)—^ ( \ 
d-2 ^4^^ \Mp) 



2p 
2p-l 



V' = 0. 

Note that these equations are consistent with the original Einstein equations ([2]) . The general solution reads 

A 



(8a) 
(8b) 

(8c) 
(8d) 



= < 



N{r) 
0(r) 
Vir) = Vo 



B 



1 - 



2p-l 

d-2p- 



Ka(-l)P2P+iC2p^, 



if b^d-3. 
if b=d-3, 



-| d-2p-l 
^ f 2p-l 



"0, 



-iV^,Cln(r) + 0o, 



if b^d-3. 



if 6 = d - 3, 



where Nao, Vq , A and <j>o are integration constants (without loss of generality, we can set A^o. 
where we have defined 



(9a) 
(9b) 
(9c) 
(9d) 

1 and (fiQ — 0), and 



B = 



2/ta(-2)^(2p- 1) / Pq \ 
(d-2)(d-2p- 1) \4ap) 



2p 
2p-l 



and C = (-2)5F^ (T'o)'""' (Aap)^^ . 



(10) 



In what follows, we consider the formalism of the grand canonical ensemble, and hence we will consider the variation 
of the action keeping fixed the temperature f3~^ and the electric potential, <I> = 4'{r+)- This variation gives bulk terms 
proportional to the constraints, some surface terms and the variation SK. As we shall see, the requirement that the 
action has an extremum, i. e. 6Ie = on-shell, will fix properly the boundary term. The condition 6Ie = implies 
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that the variation of the boundary term, which cancels out the extra surface terms coming from variation of the action 
is given by 



5K = f] 



d-2 
2k 



5K{oo) -SK{r+), 



(11) 



Since the metric solution f^{r) differs drastically for 6 7^ d — 3 and 6 = d — 3, we shall consider both cases separately. 
Interesting enough, we shall see that in both cases, the apparently divergent contributions at the infinity will cancel 
yielding to a finite and same expression in these two distinct cases. 

For 6 7^ d — 3, the variations of the fields at infinity (r —f 00) are given by 



SV\oo = SVo, 



and hence we have 



d-2 
2k 



SA 



d-2 
2k 



SB 



2p- 1 
d-2p-l 



CSVo r- 



(12) 
(13) 



(14) 



For b < d — 3, the contribution proportional to r'^~^-'' may blow up at infinity, but since the factor between round 
brackets multiplying this term identically vanishes, the variation of the boundary term at infinity yields to a finite 
expression given by 



5K{oo) = na-2SA. 
2k 



(15) 



For the case & = d — 3, which corresponds to an exponent p = [d — l)/2, we have a similar derivation. Indeed, the 
variations of the fields at infinity read 



5f\^ = -r-'''^-'^HA + Ka{-lY2P+HC^Pr-'^'^-'^hn{r), 
SV\oo = SVo, 



and hence. 



SK{oo) = (3 na- 



d-2 



— A - ((d - 2)a{-l)P2PdC^P + CSVo) ln(r) 



(16) 
(17) 



(18) 



As in the previous case, since the following equality {d — 2)a{~l)'^2'^SC^'^ + CSVo = holds for p = {d — l)/2, the 
logarithmic dependance disappears yielding to the same variation ( 15 1. We conclude that in both cases, the integration 
of the variation ( 15 1 leads to a finite expression given by 



K{oo) 



d-2 



2k 



-I3nd-2A. 



(19) 



In order to evaluate the variation of the metric function /^(r) at the horizon r_|_, we use the fact that the solution 
satisfies /^(r_|_) = together with the condition that ensure the absence of conical singularities at the horizon, that is 
= 47r//3. This condition fixes the temperature of the ensemble while the electric potential (j) and the variation 
of V are evaluated directly, 



An 

J 



(20) 
(21) 



The variation of the boundary term is easily integrated at the horizon yielding 

2tt 

K{r+) = -<^>(3Vond-2 + —nd-2ri-^- 

K 



Finally, the on-shell Euclidean action, which reduces to the boundary term K = K{oo) — K{r^), reads'^ 

Zk k 



-nd-2ri-^- 



(22) 



^ This holds up to an arbitrary additive constant that is chosen to be equal to zero by setting = in the case of flat spacetime. 
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As it will be shown below, this expression is useful in order to identify the thermodynamical quantities and in the 
study of the global stability of the black hole solutions. Indeed, the Euclidean action is related to the Gibbs free 
energy by Ie ~ f3G = f3M — /3^Q — S. This relation allows to easily identify the mass (M), the electric charge {Q) 
and the entropy (S) as 



M 

Q 
s 



f3 



OIe 
1 [OIe 

diE 

dp 



2k 



(3\d<i> 



A, 



(23a) 
(23b) 
(23c) 



As a first and direct consequence, these quantities must satisfy the first law of thermodynamics SM = TSS + ^SQ 
which is indeed the case. In a more unexpected way, these expressions (23 1 permit to derive a generalized Smarr 
formula. In order to achieve this task, we first express the event horizon radius in terms of the thermodynamical 
quantities ( 23 1 as 



2kM 



K{2p - 1) 



+ nd-2{d-2) nd-2P{d-2) 
Throughout this relation, the entropy and the temperature can be written as 



(24) 



S = At: 



2k 

nd-2[2p{d~2)Y-^ 



[2pM - (2p - 1) Q$] 



d-2 
d-3 



T 



1 /pr2d_2(rf-2) V-^ 2Mp(d-3)-2(pd-4p+l)0$ 



All \ 



{2pM -{2p- l)Q<i>)^- 



Remarkably enough, multiplying these two expressions and, after some algebraic manipulations, we obtain a formula 
for the total mass given by 



d-3 



p{d~3) 



(25) 



which is nothing but a Smarr formula. We first observe that for p = 1, the relation (25) reduces to the well-known 
Smarr formula for the d— dimensional Einstein-Maxwell black holes 0. Another interesting value of the exponent is 
given by p = 1/(4— d) for which the contribution of the charge in the formula ( 25 ) disappears and the formula is similar 
to the Smarr formula for the Schwarzschild solution. This is not surprising since in this case, the charge contribution 
in the metric (see (4b) and (|5|) becomes constant and hence, the solution resembles to the Schwarzschild metric but 
with a deficit solid angle. In fact, this critical value of p corresponds to the transition between the solutions which 
asymptote the Minkowski spacetime (type II) and those that go slowly than the Schwarzschild de-Sitter solutions, i 



fHr) 



with < — 6 < 2 (type III). This peculiar transition is also reflected in the total mass formula (25) by 



the fact that the term proportional to the charge changes its sign. More precisely, the sign of the term proportional 
to the charge in the case of solutions of type II is opposite to the sign relative of those of type III, and the transition 
is precisely operated at the critical value p = 1/(4 — d). Another value to consider is p = because of the apparent 
singularity of the formula (25). However, a detailed analysis shows that for p = the formula reduces to the standard 
Smarr formula for the Schwarzschild- (anti) de Sitter metric^ (with a cosmological constant A — Ka), 



M 



d-2 
d-3 



ST 



2anH- 



d-2 



(d-3)(d-l) + ip="^ 



where the expression of the horizon radius (24) is evaluated at p = 0. This is not surprising since at the level of the 
action, the value p — turns to be equivalent of having the Einstein action with a cosmological constant. Within this 
detailed study of the expression ( 25 ) , we have pointed out that the Smarr formula by its own reflects the different 



* The Smarr formula for AdS black holes and an extended version of the first law including variations of A have been studied in |23| . 
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FIG. 1: Mass of the black hole in terms of the event horizon radius r+ at fixed electric potential The first graph corresponds 
to the solutions with p € (|, ^^), while the second one represents the solutions with p > The last graph is identified 

with the solutions with a negative exponent p < 0. For the two first graphs (p > 1/2), the mass is a monotonous function of the 
horizon radius, with the particularity that black hole solutions with negative mass are only allowed for p > (d— l)/2. This fact 
occurs when r+ < vm, where tm corresponds to the radius for which the mass vanishes. In contrast, for p < the mass is not 
a monotonous function and has a local maximum, as it is shown in the last graph. Moreover, in the range 1/(4 — d) < p < 0, 
only black hole solutions with positive mass can be exhibited, since for this interval of p the negative mass solutions have a 
negative temperature (see Fig. 2). For the remaining range p < 1/(4 — d), the mass can be negative since its corresponding 
temperature is positive. 



asymptotic behaviors of the black hole solutions. Finally, we mention that the Smarr law can also be derived in 
two different other manners: by using the Komar integral even in the cases of the non asymptotic flat solutions, 
(see Appendix Al) or by constructing a Noether conserved current which is associated to a scaling symmetry of the 
reduced action ([t]), see Appendix A2. 

For latter convenience, we now express the mass and the temperature in terms of the horizon r+. Since we are 
working in the grand canonical ensemble, these expressions must be written in terms of the fixed state variables, 
namely the potential $ and the temperature T. A straightforward computation shows that the mass M expressed as 
a function of r+ reads 



M : 



»d-2(rf-2) 
2k 



„d-3 



d - 2p - 1 ' 



d-2p-l 



(26) 



This relation is useful to analyze the behavior of the mass with respect to the event horizon radius for the different 
classes of solutions, cf. Fig. 1. On the other hand, the square temperature is given by 



1 


dp 




47r 


dr 


r+_ 



which, after some algebraic manipulations, can be rewritten as 




D 



„2p-l 



(27) 



and where the constant D is given by 

D = -2Ka{-2)P{d - 2)-^{2p - l)^-^P{d - 2p - lfP<^>^P. 

It is important to mention that, due to the positive energy condition |14j . the constant D is positive for all the values 
of the exponent p. For clarity, in Fig. 2, we sketch the temperature as a function of the horizon radius for the different 
range of the parameter p. 



IV. LOCAL AND GLOBAL THERMODYNAMIC STABILITY 



The thermodynamic stability of a system can be considered from many different point of views depending on which 
thermodynamical variables or state functions one is considering. Usually, in order to study the stability of a system. 
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FIG. 2: Temperature against the horizon r+ for different ranges of the exponent p at fixed electric potential $. At the top left, 
the graph corresponds to the range p £ (|, ^y^), while the top right graph represents the solutions with p > The last 

graphs are respectively for p < and p £ (43^, 0). For the exponent p ranging over the set (|, ^-^), the asymptotic behavior 
of the solution is similar to the Reissner-Nordstrom black hole, and the analogy is also valid concerning the thermodynamics. 
In particular, for large value of the horizon, the temperature goes to zero while it reaches its maximum To for a finite value of 
the horizon radius vq. In addition, there exists a value rx for which the temperature vanishes, which corresponds to an extremal 
black hole. For p > [d — l)/2, the main difference with the previous case lies in the fact that there exist black holes with 
positive and negative mass, and the extremal solution is possible only for negative mass. For p < 1/(4 — d), the temperature 
behavior is quite different since it is a decreasing function of the event horizon, and as before, the extremal solution is possible 
only with a negative mass. Finally, for l/(d — 4) < p < 0, the temperature is similar to the previous case, with the exceptions 
that the mass is always positive, and the value for which the temperature vanishes does not correspond to an extremal black 
hole. 



it is common to consider small fluctuations of the state functions around the equilibrium, and since the first order 
terms vanish, the stability is only a statement about the second order variations. An equivalent manner of studying 
the local stability can be done by analyzing the sign of the heat capacity C$ at constant potential 



as well as, the sign of the electrical permittivity at constant temperature 



dQ 
9$ 



(28) 



(29) 



From these definitions, it is clear that the heat capacity (resp. the electrical permittivity) gives information about the 
thermal stability with respect to the temperature fiuctuations (resp. to the electrical fluctuations). The positivity of 
the heat capacity, C$ > 0, is a necessary condition for the system to be locally stable and, this is a direct consequence 
of the deflnition ( 28 1 and the fact that the entropy is proportional to the size of the black hole [23] . On the other 



hand, the black holes will be electrically unstable under electrical fluctuations if the electrical permittivity is negative. 
The physical interpretation of this instability is explained by the fact that in this case the potential decreases while 
the system acquires more charge, and hence the system may leave easily the equilibrium state |24| . The local stability 
of a system will be ensured only if there exists a range of the horizon radius for which the quantities C$ and ex are 
both positives. 
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FIG. 3: Specific iieat at fixed potential with p G (|, ^-j^), p > and p < 0. 



In order to express the heat capacity (28 1 as well as the electrical permittivity (29 1 in terms of the horizon radius, 
we rewrite them as 

C.-t(^Y(I!l\ ._glo._.(.-2).t'f „ ".T^-^.'-tU V PO) 

and 



d^\~^ f dQ 
dr+ Jrr. V dr+ 



However, since we have only the dependence of the charge Q in terms of the potential at constant horizon, we have 
to rewrite each factor of the above formula as 

a$ \ fdT\^^ ( dT 



Combining these last expressions, the permittivity is found to be 



p{d-2)D ^ f {2pd-6p+l)D~{2p-l){d-3)7'f-'^ 



K^^d~2p~l) " ' + y {2p-l)D-{d-^)r''^-^ 



From the above formulas (30 1 and (31), it is clear that there is a strong influence of the exponent p on the sign of 
these two expressions. Moreover, for a positive exponent p, both expressions will diverge at certain value of horizon 
radius rg, reflecting a radical change in the thermodynamical local stability of the system. Indeed, the sign of the 
speciflc heat and electrical permittivity is flipped at tq. For these reasons, it is interesting to plot C$ as well as ct 
against the event horizon r+ for the different ranges of the exponent p. 

For the range p S (1/2, (d — l)/2), there exist only black holes with positive mass. In this case, the heat capacity 
(cf. Fig. 3) for p ^ 1 have a positive and negative branches and present a vertical asymptote at tq. For black hole 
with horizon r+ G (j'Tj^'o); the heat capacity is positive, and hence the stability of the black hole is ensured for the 
thermal fluctuations. At the Reissner-Nordstrom limit, i.e. p = 1, the singularity disappears and the heat capacity 
is always negative reflecting the instability of the solution. For the electrical permittivity (cf. Fig. 4), the analysis 
is divided in two parts. For p e (1/2,1) and for the region r+ e (r^j^'o), the permittivity is positive and hence we 
conclude that the system is locally stable. In contrast, for p e (1, {d — l)/2) the system is unstable since there are 
no common regions where the heat capacity and the permittivity are both positive. For p > [d — l)/2, only small 
black hole with r_|_ G {tt,t^) and with negative mass are locally stable. Finally for p < 0, independently of the sign of 
the mass, the solution is always unstable since the heat capacity is negative and hence the system is locally unstable. 
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FIG. 4: The electrical permittivity against r+ for p G (|, 1), p G (1, ^-j^) and p > ^j^. For p < , there are two branches: 
the branch I corresponds to p > 2(s-d) while the branch II is relative to p < 2(3-d) • 

The main result is that the nonhnear electrodynamics theory considered in this paper permits the existence of locally 
stable black hole solutions. This result emphasizes the importance plays by the nonlinearity and is put in opposition 
with the local thermal instability of the standard Reissner-Nordstrom black hole. 

We now turn to the study of the global stability in order to determine whether our solutions are thermodynamically 
preferred over the Minkowski background. The Gibbs free energy G = Ie//3 is an appropriate state function to 
compare configurations in the grand canonical ensemble. For example, it is well-known that in the standard Einstein- 
Maxwell theory, the Minkowski spacetime is always favored over the Reissner-Nordstrom black hole since in this case 
the free energy of this latter is positive. In our case (cf. Fig 5) and for p < 0, since the temperature is a monotonous 
function of the event horizon radius, the free energy is a positive decreasing function and hence the Minkowski 
background is more likely than the black hole configurations. For p > (d — l)/2 and for a fixed temperature two black 
hole configurations with different free energy and size coexist. In this case, in spite of the fact that both solutions are 
less likely than the Minkowski background, there is a non-vanishing probability that the black hole with the larger 
r_(_ decays into the smaller one. In the case where the parameter p g (1/2, 1) and for the temperature T < Tq, the 
larger black holes are thermodynamically preferred over the Minkowski background, and since for T > Tq the black 
hole branch has a larger free energy than the Minkowski spacetime a first-order phase transition occurs at T = Tq- 
The situation is drastically different for p g (1, ^y^). Indeed, in this case the most probable configuration is that 
concerned with the small black hole branch, but in this case, phase transitions are not observed. 

V. CONCLUSIONS AND COMMENTS 

We have studied the thermodynamical properties of black hole arising as solutions of higher-dimensional gravity 
coupled to a nonlinear electrodynamics theory given as a power p of the Maxwell invariant. These solutions have the 
peculiarity of having different asymptotic behaviors depending mainly on the range of the exponent p, including non 
asymptotically flat spacetimes. Metrics with asymptotic relaxed fall-off has attracted much attention [25] for AdS 
gravity coupled to a scalar field with mass at or slightly above the Breitenlohner-Freedman bound [55] . This theory 
allows a large class of asymptotically AdS spacetimes where the charges can be properly defined. In the present work 
we have identified the integrations constants with the mass and the electric charge by using the Euclidean action 
in its Hamiltonian form. We have adopted the grand canonical ensemble by keeping fixed the temperature and the 
electric potential. We have shown that although the variation of the dynamical fields may diverge at the infinity, 
these divergences are canceled yielding a finite Euclidean action. From this regularized action we derive the different 
thermodynamical quantities and a generalized version of the Smarr formula is obtained. The different asymptotic 
behaviors of the solutions have been shown to be encoded by the Smarr formula. We have also proposed a derivation 
of the Smarr formula by using a Noether conserved current which results of the scale invariance of the reduced action. 
A similar derivation of the Smarr formula has been obtained in the case of scalar hairy black holes coupled minimally 
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FIG. 5: The Gibbs free energy in terms of temperature for p € (|, 1), p G (1, ^2^), p > and p < at fixed electric potential 
A first-order phase transition can be observed only in the first graph at the temperature To- 



to the three-dimensional Einstein gravity [57]. Note that in this reference the matter field vanishes asymptotically, 
but this condition is not fulfill in our case since the Maxwell potential may be divergent as r goes to infinity depending 
on the range of the exponent p. However, these divergences, which also appear in the Komar integrals, are canceled 
and a Smarr formula can be written. 

The local thermodynamic stability of the solutions have been analyzed through the heat capacity and the electrical 
permittivity. We have shown that, contrarily to the Einstein-Maxwell solution, there exist small black holes that are 
locally stable in the sense that there exists a range for the horizon radius for which the heat capacity and the electrical 
permittivity are both positive. An interesting problem to be dealt is the study of the classical dynamic stability of 
the model considered here and its possible relation with the local thermodynamic stability. 

It is worth to be noted that there exists some ranges of the exponent p for which the black hole solutions are 
preferred over the Minkowski background, and a first-order phase transition appears in the case p € (1/2, 1). This 
situation is clearly in contrast with what occurs in the standard Einstein-Maxwell case where the flat spacetime is 
always likely than the Reissner-Nordstrom solution. 



APPENDIX A: OTHER DERIVATIONS OF THE SMARR FORMULA 



1. Smarr formula from Komar integral 

Owing to the fact that the black hole solution ^ is time translation invariant, the vector field ^ — dt is a Killing 
vector field. We are going to evaluate the Komar integral for this vector over the (d — 2)-sphere at spatial infinity 
with induced metric 7, 

' rf^p.VT = - / dn^^^Vi r^-^a^a^'gu^r (Aia) 



00 



Using the Stokes theorem and the properties of the KiUing vector, the expression (Ala I can be rewritten as 



Joo Jn JT, 
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(A2) 



where S is a spacelike hypersurface covering the region between the outer horizon Ti. and the spatial infinity. Using 
the equations of motion, this last expression becomes 



(A3) 



and one recovers the Smarr formula ( 25 1 



Equating the expressions (Alb) and (A3 1, the terms proportional to the electric potential at infinity are canceled out 



Smarr formula from Noether conserved current 



The Smarr formula (25 1 can also be obtained from a Noether current density by observing that the reduced action 
([t]) , or equivalently the field equations ([8| are invariant under the following scaling transformations 



r = a r, 



— pd — 4p+l 

r{f) = a P V{r), 



—pd+Ap— 1 

(f) = o- p (j){r), 



f\f) = f\T), N{f)^a'-^N{r) 



where a is the parameter associated to the scale symmetry. Note that a similar scale symmetry has been observed 
in the case of three-dimensional scalar field minimally coupled to gravity [27j . A straightforward application of the 
Noether theorem yields the following current 



C{r) = <P 



', pd - 4n + 1 ■ 
r-p' - V 

P 



d-2 
2k 



2V 



in'Nr' 



(A4) 



which is conserved, i.e. drC{r) = 0. This conservation law can also be proved directly using the equations of motion 
(l8|. Evaluating this expression at infinity and at the horizon r_|_, one gets 



C(oo) 



d- 3 



C(r+) 



1 



'sTid-2) + P^^-^P±lQ^ 



Finally, using the fact that C is a constant, C(oo) = C(r+), the Smarr formula (25 1 is recovered 
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